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ABSTRACT 

Observational evidence shows that gravitational lensing induces an angular correla- 
tion between the distribution of galaxies and much more distant QSOs. We use weak 
gravitational lensing theory to calculate this angular correlation, updating previous 
calculations and presenting new results exploring the dependence of the correlation 
on the large-scale structure. We study the dependence of the predictions on a variety 
of cosmological models, such as cold dark matter models, mixed dark matter models, 
and models based on quintessence. We also study the dependence on the assumptions 
made about the nature of the primordial fluctuation spectrum: adiabatic, isocurvature 
and power spectra motivated by the cosmic string scenario are investigated. Special 
attention is paid to the issue of galaxy biasing, which is fully incorporated. We show 
that different mass power spectra imply distinct predictions for the angular correla- 
tion, and therefore the angular correlation provides an extra source of information 
about cosmological parameters and mechanisms of structure formation. We compare 
our results with observational data and discuss their potential uses. In particular, it 
is suggested that the observational determination of the galaxy-QSO correlation may 
be used to give an independent measurement of the mass power spectrum. 
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1 INTRODUCTION 

According to general relativity, light from distant objects 
such as quasars is affected by the large-scale distribution 
of the intervening masses, and, consequently, these act as 
gravitational lenses. A gravitational lens enlarges the solid 
angle of the source and conserves its surface brightness. So 
at the same time it brings to view faint objects (magnifica- 
tion bias) and dilutes their population density. These two 
antagonistic effects can generate either correlation or anti- 
correlation between the angular distribution of background 
and foreground objects (see e.g. Schneider, Ehlers & Faico 
1992). 

Several studies have explored the correlation between 
QSO's and galaxies at various wavelengths. For example, 
Bartelmann & Schneider (1994) found a statistically sig- 
nificant correlation between IJy sources and IRAS galax- 
ies for quasars at high redshifts. Also, Bem'tez & Martinez- 
Gonzalez (1997) reported statistically significant correla- 
tions between COSMOS/UKST galaxies and PKS radio 
quasars, and Bem'tez, Sanz & Martinez-Gonzalez (2001) 
found positive correlations using radio-loud quasars from the 
IJy and Half-Jansky samples. These results were obtained 
for small angles on the sky (of the order of 10 arcmin). Re- 



cently, Williams & Irwin (1998) investigated the correlation 
between QSO's from the LBQS catalog and galaxies from 
the APM survey and detected a significant positive angular 
correlation, ranging from scales of 10 arcmin to over a de- 
gree. For a more extended review of observational data see 
Bemtez et al. (2001), Bartelmann & Schneider (2001), and 
Norman & WiUians (2000). 

An analytical theory of the galaxy-QSO correlation due 
to weak gravitational lensing was introduced by Bartel- 
mann (1995), Dolag & Bartelmann (1997), hereafter B95 
and DB97, and Sanz, Martinez-Gonzalez & Bem'tez (1997), 
having as base the work by Kaiser (1992). We will re- 
view the theory briefiy in Section 2. The predicted corre- 
lations between galaxies and QSO's agree reasonably well 
with some observations for angles below 10 arcmin. How- 
ever, Williams & Irwin (1998) claim that for larger angular 
scales theory and observations differ by one order of mag- 
nitude for standard cosmological models such as cold dark 
matter (SCDM) and cold dark matter with a cosmological 
constant {AC DM). 

Although this conclusion has to be tested with future 
observations, it is of importance to investigate the predic- 
tions of the galaxy-QSO correlation for a large class of cos- 
mologies. The effect of bias, matter content and current ex- 
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pansion rate has to be investigated in detail. This will be the 
aim of this paper, where we keep the assumption that the 
weak lensing approximation is valid. Efltects of higher order 
terms are not investigate here. These contributions where 
estimated recently by Williams (2000), and were shown to 
increase the amplitude by less than 10 per cent for nearby 
large-scale coherent structures. Our investigations show that 
the observations of correlations between QSOs and galaxies 
on large angular scales can provide independent informa- 
tion relevant to structure formation understanding, and are 
somewhat sensitive to the values of a set of cosmological 
parameters. 

The paper is organized as follows. In Section 2 we re- 
view the theoretical formalism used. In Section 3 we present 
the cosmological models we consider in our investigation and 
study their predictions. The models include mixed dark mat- 
ter and quintessence models, as well as isocurvature and cos- 
mic string based models of structure formation. We analyze 
in detail the effect of several cosmological parameters on the 
galaxy-QSO correlation in a flat universe with cosmological 
constant and initial adiabatic fluctuations. We discuss our 
findings and some observational data in Section 4. 



2 THEORETICAL FRAMEWORK 

In this section we review the theory of the correlation be- 
tween galaxies and QSOs based on weak gravitational lens- 
ing. We follow closely Bartelmann (1995) and Dolag & 
Bartelmann (1997). However, we extend the theory by al- 
lowing for a scale- and time-dependent bias, and by using 
an improved mass power spectrum (e.g. by taking into ac- 
count the baryon density and properly evolving non-linear 
density fields). 



2.1 Angular correlation definition and calculation 

We define the angular correlation function ^gq(</>) between 
galaxies and high-redshift QSO's as 



nq 



nG{e + cj>) 



no 



(1) 



where nq and no are the QSO and galaxy densities (a bar 
over a quantity indicates its mean value), and ( ) represents 
the average over 9 and the direction of (j) (but not its modu- 
lus). Assuming weak gravitational lensing (i.e. that the mag- 
nification field contrast 5fi is small, \Sfj,\ ^ 1), and a bright- 
ness distribution of the cumulative quasar number density 
as a function of fiux S of the form nQ{> S) oc S~" , one can 
show (B95) that 



(2) 



where Sgai is the galaxy density contrast. 

The magnification fluctuations 5n derive from the in- 
homogeneities in the mass distribution, since these create 
a gravitational potential topography which deflects light, 
and therefore generates a shear field. The calculation of this 
shear field (which assumes a statistically homogeneous and 
isotropic density field and a Newtonian weak field approx- 
imation for gravity), and the results involving the power 
spectra of projections of three-dimensional fields, obtained 



by Kaiser (1992), allows one to derive the expression [equa- 
tion (12) of DB97 generalized for a scale and time dependent 
bias factor b{k, w) (Bemtez & Sanz 1999)] 
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where w is the comoving distance which here parameter- 
izes time {woo represents a redshift z = ex), and k is the 
wavenumber of the density contrast in a plane wave expan- 
sion; flm is the matter density parameter; Ps{k,'w) is the 
time evolved mass power spectrum; Jo is the zeroth-order 
Bessel function of first kind; and fxiw) is the curvature- 
dependent radial distance {— w for a fiat universe). 

The scale factor, a = 1/(1 + 2), is determined by solving 
the following integral equation for a 



w{a) — 



da' 



a'^H{a')/Ho 



(4) 



where iJ(a) = Ho [nA{l-a-'^)+nm{a-^-a-^)+a-'^Y^^ . For 
an Einstein-de Sitter Universe {flm ~ 1 and Qa = 0) the last 
equation can be easily solved, implying a(iu) = (1 — w/2)'^, 
w being measured in units of c/Ho- For a general cosmology 
equation (^ has to be solved numerically. 

The information about galaxy and QSO redshift dis- 
tribution enters through the weight functions Wg{w) and 
Gq{w), respectively. We take Wg{w) from Kaiser (1992) 
with parameters a = 1, /3 = 4 and a mean galaxy redshift 
ZG = 0.2, that is just a distribution sharply peaked at zg, 
which mimics the APM survey redshift distribution well. For 
the QSO redshift distribution we adopt a simple linear ramp 
function which starts at a high value at a minimum redshift 
zo = 0.3, and falls to zero ai z — 3.5. This also well mimics 
QSO surveys such as LBQS and PKS. From this distribu- 
tion the QSO weight function Gq(w) is calculated following 
DB97. Note however that the exact shape of these distri- 
butions is not crucial for our results, and its role is already 
discussed in B95. 



2.2 Mass power spectrum 

We define the linear mass power spectrum as usual 

A\k,^)^^Psik,n.)^A(^) 

where A is the normalization factoiQ ; n is the initial power 
spectrum index (for a Harrison-Zel'dovich spectrum n = 1); 
T{k, w) is the transfer function; and D{w) is the growth 
function which we define as 



T\k,w)D^(w), (5) 



(6) 



* A = (5 if we normalize to COBE (Bunn & White 1997). An- 
other possibility is to normalize to the cluster abundance (Viana 
& Liddle 1999), which would be in principle the most appropriate 
for our calculations, because it normalizes the power spectrum on 
scales that are relevant to weak lensing, while a normalization to 
the cosmic microwave background anisotropics does the normal- 
ization on much larger scales. 
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where a — a{w), and g{a) is the linear growth suppres- 
sion factor which is well approximated by (Carroll, Press 
& Turner 1992; Lahav et al. 1991) 



5(a) 



2 

+ (1 

^71 



^/'(a)-f7A(a) 
1 



J^m(a) = — ( 7f 



2 



l + ^r!A(a; 



Ha (a) = riA 



(7) 
(8) 
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For an Einstein-de Sitter Universe g{a) — 1 (no suppression 
of gravitational clustering) , and the growth function reduces 
to the scale factor, D{w) — a{w). 

For CDM we use the transfer functions given by 
Bardeen et al. (1986), calculated with the shape parame- 
ter by Sugiyama (1995), T = Q,mhexp[—^lbii + V2h/^m,)], 
where Qb is the baryon energy density. 

A reahstic model for the mass power spectrum has to 
consider non-linear effects in the evolution of the density 
fields, important at small scales. One way to do this is by 
computer- intensive simulations, another is to use an approx- 
imate analytical form such as a scaling equation mapping the 
linear spectra to the non-linear one (Hamilton et al. 1991, 
Scranton & Dodelson 2000) 



A%L{kNL,'w) = fNL (Ai(fci,,w)) 



1/3 



(10) 



fu) 



Expressions for the function fNL(x) are given by Peacock 
& Dodds (1996) and Ma (1998), hereafter PD96 and Ma98. 
In the linear regime [x <^ 1) fNL{x) ~ x, during the stable 
clustering regime {x ^ 1) fNL{x) oc x^^'^, and in the inter- 
mediate region (x « 1) fML{x) is extracted from N-body 
simulations. 



2.3 Galaxy biasing 

The simplest bias relation one can assume is linear and de- 
terministic, so that one could define the bias as 



b = 



where 



(i?) = / W^{kR)A\k) 



dk 
k 



(12) 



(13) 



and W{x) = 3 [sin(x) — a; cos(2;)] /a;^. Evaluated at i? = 
8h~^ Mpc, this gives erg. The subscripts 'gal' and 'mass' refer 
to the galaxy and mass power spectra. 

However, the linear deterministic bias is a too simplis- 
tic model (Dekel & Lahav 1999). The most general biasing 
between galaxy and mass density contrasts that one can 
imagine is a non-linear stochastic bias which can be time- 
and scale-dependent: 



Sgat = b{S)S + €, 

where e is a random field. 

From equation (^) we get that 



(14) 



= {5t,{e)b[S{e + $)]S{9 + $))+{s^^{e)eie + $)) .(15) 

If we make the reasonable assumptions that e is independent 
of Sjj. and has zero mean, then the last term on the right is 
clearly zero^ , so a stochastic component of the biasing does 
not alter the expected value for the galaxy-QSO correlation. 

To be most general, we further consider a scale- and 
time-dependent bias factor, which we assume can be written 
as b{k, w) = bs{k)bt{w). The scale-dependent part of the bias 
is defined as 



bs{k) 



Ps{k) 



(16) 



where Pgai (k) is an empirical galaxy power spectrum. We use 
(unless indicated) the spectrum obtained from the APM sur- 
vey (Gaztanaga & Baugh 1998). The time-dependent part 
of the bias can be modeled as (Tegmark & Peebles 1998) 



bt{w) = 



^(1 - DY - 2(1 - D)robo + bg 
bo 15 



(17) 



where bo is the linear bias factor today (equation |12| ) , ro is 
the dimensionless correlation coefficient between the distri- 
butions of mass and galaxies, and D — D{w) (equation |^). 
However, we found that the time dependence is negligible 
(even considering a galaxy population at a higher redshift 
and a strong time dependence of the bias the effect is very 
small), so for simplicity one can neglect it. 



3 RESULTS 

In this section (except for Section 3.4) we explore the de- 
pendence of the galaxy-QSO angular correlation on several 
cosmological parameters, assuming structure formation from 
initial adiabatic fluctuations. 



3.1 Three historical models 

Fig. |l| shows our results for three cosmological models that 
were in fashion at some point in history, the 'standard' cold 
dark matter {SCDM: Q,m — 1, h — 0.5), a flat universe 
with a cosmological constant {ACDM: flm = 0.3, Ha ~ 0.7, 
h — 0.7), and an open universe (OCDM: fim = 0.3, SIa = 0, 
h — 0.7). It uses the COBE normalized mass power spec- 
trum for these models, the necessary bias to match the APM 
galaxy power spectrum (equation jl6j), and the calculated 
galaxy-QSO angular correlation (equation [^). The three 
cosmological models are clearly differentiated by the pre- 
dicted ^GQ{<t>)- 

A normalization of the power spectrum to the clus- 
ter abundance simply rescales our results for ^gq{4') by 
{o'8)ciuster/{o'8)coBE- The resulting correlations decrease 
for SCDM, increase for OCDM, and remain the same for 
ACDM, making the differentiation among the models less 
evident. 

The latest CMB observation (Bernardis et al. 2000, and 

t While this is certainly true for averages over large angular re- 
gions, for patches smaller than the typical scales of variation of 
the fields involved this term could contribute because of statistical 
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Hanany et al. 2000) point to a flat universe, as favored by 
most inflationary models. The result of these observations, 
and the suggestion of the existence of a dark energy com- 
ponent by supernova observations, strongly induce one to 
take n = Qm + = 1- This flat universe with low matter 
density (most of it coming from some sort of cold dark mat- 
ter), dark energy (cosmological constant or quintessence), 
and adiabatic fluctuations with an initial spectrum nearly 
scale-free, constitutes the most popular cosmological model 
today. We will concentrate on this model, searching for the 
efi^ects of variations of the cosmological parameters around 
the most likely values. 

3.2 Adiabatic flat universe with A 

The matter density parameter flm enters the calculation of 
^gq{<I>) at several levels. It is a multiplicative factor of the 
whole expression in equation (^, it determines the general 
geometric scaling of the problem through the scale factor, 
and it also enters in the mass power spectrum (and conse- 
quently in the bias, too) through the shape parameter, the 
growth function, and normalization. 

We investigate the dependence of the galaxy-QSO cor- 
relation on the matter content of a flat universe in Fig. ^. As 
would be expected, a higher matter content generates more 
power of the mass spectrum on small scales, which results 
in a higher correlation between galaxies and QSO's (mainly 
at small angles). 

Recent experiments give extra support to the possibil- 
ity that neutrinos are massive. This necessarily implies a 
hot dark matter (HDM) component to the matter content 
of the universe, because theory predicts a cosmic neutrino 
background of temperature 1.95 K. A pure HDM universe is 
ruled out, because it would not have the observed amount 
of structure on small scales, as they are suppressed by free 
streaming. However, a mixed dark matter scenario (MDM) 
with a small fraction of HDM is possible. 

We calculate ^gq (0) in a MDM universe using the linear 
mass power spectrum derived by Ma (1996) and obtain the 
full non-linear time evolved spectrum using the prescription 
by Ma (1998), which properly takes into account the time 
evolution of the hot and cold dark matter components (the 
PD96 prescription does not give the correct result in this 
case). Fig. ^ shows our results. The larger the fraction of 
HDM over CDM, the larger the suppression of structures on 
small scales. A large neutrino energy density fli, results in 
less power at large k. The consequence on ^gq(0) is a lower 
correlation at small angles (and a flatter general slope). 

The baryon density of the universe is fairly well deter- 
mined (Buries et al. 1999) and represents a small fraction 
of the total matter. In our calculations Qb enters through 
the shape parameter F of the transfer function. This is just 
an approximation, because it does not account for baryonic 
oscillations (Eisenstein & Hu 1998), but it is good enough 
for low values of Qt- We studied the influence of a baryonic 
content varying from to 8 per centof the critical density 
on the galaxy-QSO correlation (see Fig. A higher baryon 
density implies less power on small scales for the mass spec- 
trum and a lower value of £,GQ{<t>) on small angles. 

Inflation favors an initial power index n « 1 for the 
spectrum of density fluctuations, but some degree of 'tilting' 
is possible depending of the inflationary theory. We look at 



the effect of modifying the initial power spectrum index by 
simple tilting at Fig. ^. A higher n implies a steeper ^gq{4>), 
since it induces more power on small scales. 

The effect of a Hubble parameter in the range 0.5 < h < 
0.8 on the predicted S,gq{(I>) was also considered. A higher 
value for h implies a higher correlation, as also previously 
shown in DB97. 

3.3 Quintessence 

One alternative for the cosmological constant as dark en- 
ergy is a dynamical and spatially homogeneous (on cluster 
scales) field, or, as it has been called, quintessence (see e.g. 
Wetterich 1988, Ratra and Peebles 1988). Such a field has 
negative pressure and an equation of state p — ujqp with 
coefficient in the range —1 < ujq < 0. We implement a 
quintessence model (QC'DM) with constant ujq by general- 
izing the time dependent Hubble expansion rate to 

H{a) = Ho [f2Q(a-^("«+')-a-')+n™(a-='-a-')+a'Y''' ,(18) 

and by using the results of Ma et al. (1999) for the con- 
struction of the fully evolved non-linear mass power spec- 
trum. Note that loq = — 1 reduces everything to the case of 
a cosmological constant. 

Our results are shown in Fig. ^ A greater cuq (closer to 
0) implies a lower amount of large-scale structures over small 
ones in the relevant range (at much larger scales the clus- 
tering of the quintessence field becomes important), which 
results in a steeper slope for £,gq{4')- The same conclusion 
is valid if we normalize the QCDM spectra to the cluster 
abundance, which also depends on ujq according to Wang 
& Steinhardt (1998). Note that, because the larger scale 
structures are the most affected by the coefficient of the 
quintessence equation of state, the effect on the correlation 
is more dramatic at large angles. 

3.4 Alternatives for structure formation 

We also investigated alternatives to adiabatic fiuctuations 
as mechanism for structure formation: isocurvature fluctua- 
tions and the cosmic strings scenario. 

A general treatment of initial perturbations allows the 
existence of isocurvature modes (Bucher, Moodley, & Turok 
2000). We consider an isocurvature CDM model for the mass 
power spectrum using the transfer function by Bardeen et 
al. (1986) and the spectral index obtained by Peebles (1999). 
In this case we use PD96 to obtain the non-linear spectrum 
and normalize it to the cluster abundance. 

Except for inflationary models, topological defect mod- 
els are the only known way of seeding structure formation 
(see e.g. Vilenkin and Shellard 1994, Hindmarsh and Kibble 
1995, and Brandenberger 1994 for recent reviews). We ex- 
plore this scenario for structure formation through the mass 
power spectrum generated by a network of cosmic strings 
(Avelino, Caldwell & Martins 1997). Note that the formal- 
ism used is insensitive to deviations from Gaussianity. Such 
deviations are a common feature of cosmic strings scenarios. 
Again, the non-linear spectrum is obtained by PD96 pre- 
scription and is CMB-normalized, which implies ag, — 0.46. 
A normalization to the cluster abundance requires erg = 0.81 
[see Avelino, Wu & Shellard (2000); a somewhat smaller nor- 
malization was obtained by van de Bruck (1999)]. 



© 0000 RAS, MNRAS 000, 000-000 



Galaxy- QSO correlation 5 



Both scenarios are in serious trouble if considered as the 
only mechanism responsible for structure formation. How- 
ever, hybrid models are still allowed and maybe even favored 
by the first results from the Boomerang data (Bernardis et 
al. 2000) which yields an indication of a low amplitude for 
the secondary acoustic peak. See Bouchet et al. (2000) , Bat- 
tye & Weller (2000), and Avelino, Caldwell & Martins (1999) 
for considerations on hybrid models with topological defects, 
and Enqvist, Kurki-Suonio & Valiviita (2000) for limits on 
isocurvature fiuctuations from CMB anisotropies. Our re- 
sults for these two alternative models for structure forma- 
tion are compared with the adiabatic case in Fig. ^. Both 
predict a steeper slope for £,gq{4i) than adiabatic fluctua- 
tions. This is because their power spectra have more power 
on smaller scales relative to large ones than the adiabatic 
spectrum. 



4 DISCUSSION 

We studied in detail the role of several cosmological param- 
eters and models of structure formation on the prediction of 
the angular correlation between galaxies and QSOs induced 
by weak gravitational lensing. 

In a log- log plot ^GQ (0) can be roughly characterized by 
an amplitude and a slope. The amplitude is closely related 
to the normalization of the mass power spectrum, and the 
slope to its shape. More power on small scales (large k) im- 
plies a steeper slope. This is a reflection of the fact that the 
amount of small scale structure is responsible for the galaxy- 
QSO correlation on small angles and larger structures are 
responsible for the correlation on larger angles. This sim- 
plicity of behavior results in a relatively large degeneracy 
of 5 GQ (<^) as a function of the cosmological parameters and 
models analyzed. In brief, a higher value of ilm, n, or h has 
nearly the same effect on ^gq('^) as a lower value of fl^, 
or Qf This effect is a higher correlation on small angles. 
Roughly, this represents a steeper slope for ^gq(<^), which 
can be also achieved by shifting luq closer to zero, or by 
introducing an isocurvature or cosmic string component. 

The main lesson given by the use of different theoretical 
mass power spectra is that the galaxy-QSO correlation at 
different angles contains information about the amount of 
structure at different scales. 

Our calculations fully incorporate galaxy biasing. This 
requires the use of an empirical galaxy power spectrum. We 
choose to use the spectrum obtained from the APM sur- 
vey, but this choice is not unique. Spectra from other sur- 
veys are available and could be used as well. These spectra 
do not all agree (see Einasto et al. 1999), allowing some 
freedom for the galaxy-QSO correlation prediction, because 
our ^gq(<^) calculation can be seen as an integral over an 
effective mean spectrum -^Z Pgai{k)Ps{k). In a similar way, 
the use of a power spectrum of cluster of galaxies allows 
for a prediction of the correlation between galaxy clusters 
and QSOs. To illustrate this point we use the spectrum ob- 
tained from the Abell-ACO clusters of galaxies (Miller & Ba- 
tuski 2001) as Pgai{k). We show the effect of this alternative 
choice at Fig. ^ which gives a higher amplitude for ^gq (</>), 
roughly given by a factor {as)AbeU-ACo/{^8)APM ~ 4, 
where {(ys)Abeii-ACO ~ 3.2 is the rms fluctuation for the 
Abell-ACO spectrum, and {a8)APM = 0.82. The correlation 



curve found using the Abell-ACO spectrum is also slightly 
steeper than that using APM, because the first spectrum is 
fiatter than the second at large k [PAbeU-ACo{k) oc k~^'^ 
and PAPM{k) oc fc~^'^]. 

We also plot in Fig. ^ some points obtained from obser- 
vations. A direct comparison among the different data sets is 
not possible, because they were obtained from different QSO 
and galaxy populations. For the same reason a strict com- 
parison between these observational points and the theoret- 
ical predictions shown is not possible. That would require 
an individual calculation for each of the data sets, which 
would take in account the different redshift distributions for 
QSOs and galaxies, and also the particular power spectrum 
of the foreground population (galaxies or galaxy groups). 
This would be necessary because when selection criteria are 
applied to a foreground object catalog the effective power 
spectrum of the resulting subset is possibly different from 
the spectrum of the whole population. A rigorous analysis 
of the data would also require considerations about observa- 
tional particularities such as the possible effects of dust and 
systematics. That is not our aim here. Nevertheless, a qual- 
itative agreement in behavior between data and predictions 
can be glimpsed. 

Another way to analyze the results from the several 
groups (avoiding the messiness of noise data) is to capture 
the essence of each data set by adjusting a power law to 
each one. We find the best fit of Cgq(0)/(s - 1) = A(^l')^ 
to the data by the least squares method (see Table llj) . To 
compare with the predictions shown in Fig. ^ we can take 
the amplitude A as being the value of ^GQ{(t>)/{s — 1) at 
= 1' (for the Abell-ACO result A ^ 0.04 and for APM 
A « 0.008), and power index B as the average bi-log slope 
of Cgq(0)/(s — 1), which gives B ~ —0.8. Note however 
that for (j) < 1' the weak lensing assumption is not satisfied 
(strong lensing becomes important), so that the predictions 
for £,gq{4>) are underestimated in this range. 

Bem'tez & Martinez-Gonzalez (1997) (BMG97) found 
the galaxy-QSO correlation for angles up to 15 arcmin us- 
ing radio loud QSOs from the PKS catalog and also using 
optically selected QSOs from the LBQS catalog. For this 
last group there is a large dispersion of data and points 
with negative correlation on small angles were found; the 
overall result is a null correlation (dust and selection bias 
were mentioned as possible explanations). The large error 
bars obtained for the two quasar samples make these obser- 
vations compatible with predictions. Benftez et al. (2001) 
(BSMGOO) obtained a somewhat higher correlation on the 
same anglular scales using quasars from the Half-Jansk and 
IJy samples. 

Williams & Irwin (1998) (WI98) also used a subsample 
of the LBQS catalog to determine the galaxy-QSO correla- 
tion, but for larger angles (7 to 300 arcmin). For angles up 
to nearly 100 arcmin they find a much higher correlation 
than the predicted one, and for angles larger than that the 
correlation is dispersed around zero, being compatible with 
predictions. Croom & Shanks (1999) (CS99) reinterpreted 
the anti-correlation between faint QSOs and galaxy groups 
found by Boyle, Fong & Shanks (1988) as being due to grav- 
itational lensing. A flat QSO number count (s = 0.78) was 
used to explain the negative correlation, but a surprisingly 
high amplitude is obtained when ^gq{4')/{s — 1) is calcu- 
lated. This high amplitude would be in concordance with a 
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higher normalization of the power spectrum of the galaxy 
groups. Other works also seek to find the observed galaxy- 
QSO correlation, but the general conclusion that one can 
reach is that up to now (,GQ{(t>) is an observationally not 
very well determined quantity. 

The use of much larger catalogs (e.g. the Sloan Digital 
Sky Survey with « 10^ QSOs, and the 2dF survey) may 
allow a more precise and accurate measurement of ^gq(0) 
which in turn may constitute an independent way for the de- 
termination of the mass power spectrum, and a valuable tool 
for investigating structure formation and large-scale struc- 
ture. 
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Figure 1. Three cosmological models. Non-Linear mass power 
spectrum (top graph), bias (central graph), and angular correla- 
tion Cgq('^)/(*~1) (bottom graph). The internal plot is the result 
for the mass power spectrum normalized to the cluster abundance. 
Solid lines are for SCDM, Cl^ = 1, h = 0.5 {as = 1.1); dashed 
ones for ACDM, dm = 0.3, Qa = 0.7, h = 0.7 {as = 1.0); and 
dotted ones for OCDM, Clm = 0-3, Qa = 0, h = 0.7 {as = 0.46). 
Other parameters are Clf, = 0.019//i^, n = 1, non-linear approxi- 
mation by PD96. 



Figure 2. Dependence on matter density in a flat universe with 
cosmological constant (Qrn+^A = 1)- Solid lines are for Qm = 0.6 
{as = 1-4); dashed ones for = 0.4 (erg = 1.2); dotted ones 
for Urn = 0.3 {as = 1.0); and dot-dashed ones for Um = 0.2 
{as = 0.72). Other parameters are Q(, = O.OlQ/ft^, h = 0.7, n = 1, 
non-linear approximation by Ma98. 
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Figure 3. Dependence on hot dark matter density. Solid lines 
are for = {as = 1-0); dashed ones for = 0.1 {erg = 0.84); 
and dotted ones for = 0.2 (erg = 0.81). Other parameters arc 
{arn,^iA,^b) = (0.3,0.7,0.019//i2), h = 0.7, n = 1, non-linear 
approximation by Ma98. 



Figure 4. Dependence on baryon density. Solid lines are for Qij = 
(o"8 = 1.3); dashed ones for Qf, = 0.04 (erg = 1.0); and dotted 
ones for = 0.08 (eg = 0.79). Other parameters are Qm = 0.3, 
Ha = 0.7, h = 0.7, n = 1, non-linear approximation by Ma98. 
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Figure 5. Dependence on the primordial spectral index. Solid 
lines are for n = 1.15 (as = 1-4): dashed ones for n = I (eg = 1.0); 
and dotted ones for n = 0.85 (trg = 0.75). Other parameters are 
{arn,^iA,^b) = (0.3,0.7,0.019//i2), h = 0.7, non-linear approxi- 
mation by PD96. 



Figure 6. Quintessence. Solid lines are for ujq = — 1 {as = 1.0); 
dashed ones for ujq = —0.7 (erg = 0.9); dotted ones for luq = —0.5 
(erg = 0.76); and dot-dashed ones for luq = —0.3 (erg = 0.54). 
Other parameters are {nm,iiQ,^b) = (0.3, 0.7, 0.019//i2), h = 
0.7, n = 1, non-lineax approximation by Ma at al. (1999). 
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Figure 7. Alternative models for structure formation. Solid lines 
are for adiabatic fluctuations (n = 1, erg = 1.0); dashed ones for 
isocurvature fluctuations (n = —1.8, erg = 0.99); and dotted ones 
for cosmic strings {as = 0.46). The thin dotted line corresponds 
to the normalization of the cosmic strings spectrum to the cluster 
abundance (ug = 0.81). Other parameters are Clm = 0.3, Qa = 
0.7, = 0.019//j^ (cosmic strings results do not incorporate 
baryons), h = 0.7, non-linear approximation by PD96. 




Figure 8. Predictions and some observational points. Solid lines 
are the theoretical predictions for a fiat universe with cosmological 
constant and adiabatic fluctuations (Clm = 0.3, Ha = 0.7, f2(, = 
0.019//i^, n = 1, h = 0.7, as = 1.0). The upper curve uses the 
power spectrum from Abell- ACO (galaxy clusters) , and the lower 
curve the APM survey (galaxies). Observational points arc also 
plotted as an illustration. Some error bars would explode in the 
lower part if plotted, so we do not show them. The points close 
to the lower axis represent values of ioQW/i^ ~ 1) equal or less 
than zero. Note the caveats listed in the main text concerning the 
comparison among different data sets and theoretical predictions. 
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Table 1. Best fit to S.GQ^'t') / ~ ^) — ^(</'/l')^i with B < 0. is the angular range covered in arcmin, and A'^ is the number of points. 
The uncertainties are given by the projection of the la confidence level contour, and (...) represents a large value outside the investigated 
parameter space. Values for s were taken from the respective papers. 
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